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Abstract. Within the Lindblad formalism we consider an interacting spin chain coupled locally to heat 
baths. We investigate the dependence of the energy transport on the type of interaction in the system as 
well as on the overall interaction strength. For a large class of couplings we find a normal heat conduction 
and confirm Fourier's Law. In a fully quantum mechanical approach linear transport behavior appears to 
be generic even for small quantum systems. 

PACS. 05.60.Gg Quantum transport - 05. 30.-d Quantum statistical mechanics - 05.70. Ln Nonequilibrium 
and irreversible thermodynamics 



1 Introduction 

Heat conduction in condensed matter, especially in in- 
sulators, is a long standing problem which has attracted 
renewed interest recently. In the classical understanding, 
heat conduction in insulators results from phonon scat- 
tering. This process is an anharmonic effect of the lattice. 
There are two different scattering types, the normal or re- 
processes and so called U-processes. N-processes are mo- 
mentum conserving and cannot affect the heat conduction 
■ of the solid state. Only the second type of scattering, the 
U-processes, should give rise to a finite heat conductivity: 
In such processes momentum is conserved only modulo 
a reciprocal lattice vector |15l9ll(Jl7) . For very low tem- 
peratures these U-processes rapidly die out and the heat 
conductivity diverges. In this case only impurities may 
limit the heat conductivity. 

Recently, a series of articles have addressed heat con- 
duction and Fourier's law in one dimensional systems. It 
has been found that regular heat conduction (non vanish- 
ing local gradient of temperature) does not occur naturally 
in the Hamilton models considered: Most of these have 
been classical one-dimensional systems coupled via differ- 
ent interactions |5I1I14| . Note that because of the lack of 
phonon scattering, a harmonically coupled chain cannot 
have a finite heat conductance ^21- It is necessary to take 
an anharmonic potential into consideration. 

Heat conduction is a problem of non-equilibrium ther- 
modynamics in which equilibrium properties hold at most 
locally. The equilibrium thermodynamics of small systems 
can be given a quantum mechanical foundation 4 2 3 and 
therefore it is tempting to expect non-equilibrium behav- 
ior also to follow directly from quantum mechanics. In- 
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deed, in one-dimensional quantum systems strong indi- 
cations have been found for a normal heat conduction 
|12lllj . But it is not clear yet which conditions a system 
has to fulfill in order to exhibit local thermodynamical be- 
havior. Especially the dependency of heat conduction on 
the type of internal couplings within the system have not 
been investigated in full detail. There are systems show- 
ing a normal heat conductance |1 211 f| and others where 
the conductivity diverges [H] (vanishing temperature gra- 
dient). Unfortunately, even the concept of local thermo- 
dynamical quantities such as local temperature or local 
energy easily becomes ambiguous. Nevertheless it is in- 
evitable to discuss these questions carefully in the context 
of heat conduction models, since the results depend cru- 
cially on the underlying concepts. 



To analyze the conditions of thermodynamical behav- 
ior in small systems we study a one-dimensional quantum 
system coupled to heat reservoirs. The building blocks 
(here N = 4 — 6 spins) of the system are internally coupled 
with different types of interactions. Traditionally, thermo- 
dynamical behavior is expected for large systems (N — » 
oo), smaller quantum systems are thought to tend to de- 
viate from thermodynamical behavior. However we could 
find normal heat conduction even for very small systems 
(N = 4). We consider a Forster coupling alone and to- 
gether with a non-resonant coupling and we introduce a 
random coupling, too, to test the generality of thermody- 
namical behavior. Additionally, we investigate the trans- 
port behavior with respect to different coupling strengths. 
For these various cases we find diverse transport behavior 
and therefore means to specify conditions a system has to 
fulfill in order to show a normal transport behavior. 
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2 Model 

We consider an open quantum system, a chain of N sub- 
systems with n levels each, with nearest neighbor inter- 
actions (cf. 0ED- The hamiltonian H of the system thus 
reads 



A' 



N-l 
0=1 



(1) 



Here, the terms of the first sum represent equidistant n- 
level subsystems. The second term accounts for pair-inter- 
actions between two adjacent subsystems, with an overall 
coupling strength A. For A to characterize the absolute 
interaction strength it is necessary to normalize the dif- 
ferent interaction types by /, the mean absolute value of 
interaction matrix elements, i. e. 



N-l 



I 2 



(2) 



For n = 2 (spins) it is possible to use as a set of basis 
operators the Pauli spin operators &i (i — Q,x,y,z). In 
terms of these operators the local hamiltonian of a sub- 
system /j, with an energy spacing AE — 1 can be written 

as 

H^\p) = \cr z {p)- (3) 

The spins are coupled with alternative types of next neigh- 
bor interactions (see [Hj): a non-resonant diagonal inter- 
action 



H 



(NR.) 



i = Cnr o-z(fJ-) ® 0"«(M + 1), 



(4) 



a resonant energy transfer interaction (Forster-Coupling) 



H 



S+i = -y {°x (m) ® + 1) 

+ <7,,(/x) ®<7 2/ (/X+ 1)), (5) 



where Cnr and Cf can be used to adjust the relative 
strength, or a totally random next neighbor interaction 



n 



(R) 



(6) 



with normal distributed random numbers p^- (variance 1). 
Note that Pij is independent of p (no disorder). The ran- 
dom interaction is supposed to model " typical interac- 
tions " without any bias. 

To model the influence of a heat bath, we use the Lind- 
blad formalism: The standard Liouville equation is supple- 
mented by an incoherent damping term (see, e. g. , [H]): 



dp 
dt 



-i [n,p] + c^ t 



= Cp. 



(7) 



For the local coupling of spin p of the chain to a bath we 
expand in terms of raising and lowering operators: 



Wi- 



tt', 



P<7+ - pcr+a^ - &+&^ p) 



0->l / -H 



(8) 
10) 



where the first term describes a decay from |1) 
with rate Wi_»o an d the second from |0) — > |1) with 
Wo->i < Wl_,.o- The properties of the environment (bath) 
only enter via these two rates. 

If the bath was brought in contact with a single spin, 
the spin would relax to the stationary state: 



1 



Pstat 



Wi- 



w^ a 

WW 



(9) 



Interpreting this state as the result of thermal equilibrium, 
where T is defined via the Boltzmann distribution 



W . 



W 



l->0 



(10) 



Operationally, we can associate the temperature T with 
the temperature of the bath T(B) = T measured in units 
of AE. With fTr))l and the normalization W\-, + Wo-»i = 
\( B ^ it is possible to rewrite the rates in terms of the cou- 
pling strength A' B ) and the bath temperature T(B): 



WW = ^A (B) , W^ 

1 + c T < B > 



1 + c f < B > 



(11) 



To measure the local temperature of a spin p in the 
chain we use the mean excitation energy 



T(p) = Tr{p(p)H iioc \p)}. 



(12) 



of the respective spin (0 < T(p) < 0.5). H {loc \p) and thus 
T(p) are defined in units of AE, p(p) is the correspond- 
ing reduced density operator. This T{p) is a well-defined 
local quantity, irrespective of further couplings |16) . For 
consistency reasons we define the temperature of the bath 
now as the mean energy a single spin coupled to this bath 
would have 



T(B) 



Tr{p(l)W( loc )(l)} 
Wo-,1 _ 1 

X(B) „ — 1 — 

^ l_|_ e T(B) 



(13) 



Rather than applying iterative numerical methods (e.g. 
Runge-Kutta) to integrate the Liouville equation, we diag- 
onalize the Liouville operator C (see eq. 0) of the whole 
open system. We thus obtain the exact solution of eq. Q 
at all times and without iterative numerical errors. 



3 Single bath 

In order to investigate whether full thcrmodynamical equi- 
librium may be reached for a chain model coupled only 
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Fig. 1. Temperature profile: N=4 spins coupled by a random 
interaction (of internal coupling strength A = O.Of , 0.3, 10); 
first spin fi = 1 is coupled to a bath with coupling strength 
A (B) = 0.01 and a temperature T(B) » 0.27 (line in the box 
marks the temperature of the bath). 



at one border to a heat reservoir as mentioned in the 
last section, we consider a system of N = 4 spins, where 
the first spin (p = 1) is coupled to a thermal bath. The 
temperature of the bath is still defined by the two con- 
stants W^i.+o and Wo— >i as mentioned in the last section 
(T(B) = 1 =>• T(B) » 0.27 cf. l(T5)ll. We choose the cou- 
pling strength of the first spin to the bath to be weak 
(A^ = 0.01). 

Firstly the coupling between nearest neighbor spins is 
taken as the random interaction Ti.( R \ From the station- 
ary state of the system we compute the reduced density 
operators of each spin p(p) and the appropriate tempera- 
tures based on eq. i|12|l . The final state is independent of 
the initial state but not of the internal coupling strength 
A. For very small A = 0.01 the temperature profile is flat 
(see FIG. ^1 while for A = 0.3 the spins do no longer have 
the same temperature. In case of very strong interactions, 
A = 10, the individual spins are in a totally mixed state 
and the whole system is highly entangled. 

As a locally defined thermal equilibrium state for such 
a system one may expect the (separable) product state: 



Peq = /Ostat (/■»)• 



(14) 



To measure the distance between the actual final state pf 
and this expected state p eq we use 



A = Tr{(p cq -p f ) 2 }. 



(15) 



In FIG. we vary A = 0.01. ..1.3 for different random in- 
teractions and evaluate the corresponding distance to the 
equilibrium state p eq in a N = 4 spin chain. For decreas- 
ing A we find that the resulting pf is indeed coming closer 
to the expected stationary state. For very strong inter- 
actions we find that the final state is very close to the 
totally mixed state (D 1 = Tr{(/5 oq - ^vi) 2 } = 0.073). 
Qualitatively all the random interactions show the same 
behavior. 

However, if we considered instead of the random inter- 
action the case Cnr = Cf = 1 (Forster and non-resonant 
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Fig. 2. Distance Df between the expected equilibrium state 
/3 eq and the final state pf in dependence of interaction strength 
A for six different random interaction in a N — 4 spin chain. 



coupling) we would not get any dependence on the in- 
ternal interaction strength A: Irrespective of A all spins of 
the system have exactly the same temperature as the heat 
bath. 

In conclusion we always get thermodynamical behavior 
for the special case of Forster and non-resonant coupling. 
Generally (random interaction), thermodynamical behav- 
ior can be expected in the weak coupling limit. 



4 Two Baths 

Now we go on to study non-equilibrium properties. For 
this purpose we consider again the chain of N = 4 spins, 
but now coupled at both ends to separate baths of differ- 
ent temperature but the same A^ B ) = A^ B1 ) = A^ B2 ^. For 
the higher temperature bath Bl we set Wo_>i — 0.4A^ B \ 
VFi^o = 0.6A < - B - ) , which corresponds to a bath tempera- 
ture in our definition of T(B1) = 0.4. This bath is coupled 
to the first spin (p = 1) in the chain. The lower tem- 
perature bath B2 (T(B2) = 0.2 with W ^i = 0.2A( B >, 
Wi^q — 0.8A( B )) is coupled to the last spin (p = 4) in 
the chain. Again we assume a weak coupling, A^ B ) = 0.01, 
of the spin system to both baths and vary the internal 
coupling strength of the system A. 

Here we investigate all three different types of inter- 
nal couplings between the spins to test the dependence 
of the interaction. Firstly we consider the coupling due to 
an energy transfer only (Forster interaction Ti.^ F \ Cf = 1, 
Cnr = 0). For a strong coupling within the chain, A = 1, 
all spins get the same averaged temperature between those 
of the heat baths, independently of the initial state. In 
case of a weak internal interaction strength, A = 0.01, 
the two borderline spins are drawn to the temperatures of 
their bath (T(l) = 0.302 and T(4) = 0.298, respectively 
cf. FIG. |2l • The two spins in the middle are exactly on 
the same average temperature T(2) = T(3) = 0.3 inde- 
pendently of the interaction strength (see FIG. EJ. Thus, 
in the middle of the system the conductivity diverges and 
Fourier's Law is not fulfilled. This is a well known behavior 
for systems with only energy transfer coupling |5I13| . 

Now we couple the subsystems additionally with a non 
resonant interaction, which in itself does not give rise to 
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Fig. 3. Temperature profile: spin chain (N — 4) with only 
internal next neighbor Forster coupling (Cnr = 0, Cf = 1) 
coupled to heat baths at both ends (lines in the bath boxes 
mark the appropriate temperature of the bath T(B1)=0.4, 
T(B2)=0.2); internal and external coupling strength A = 
A< B » = 0.01. 



Fig. 5. Temperature profile: open system with internal ran- 
dom next neighbor coupling for increasing internal coupling 
strength A (A (B) = 0.01); bath temperatures T(B1)=0.4, 
T(B2)=0.2. 
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Fig. 6. Distance Df between the expected equilibrium state p oq 
and the final state pi as a function of the interaction strength 
A for two different random interaction in a N = 4 spin system. 



Fig. 4. Temperature profile: open system with N = 4 and 
N = 6 spins coupled with internal next neighbor Forster and 
non resonant coupling (Cnr = f, Cf = f); internal and ex- 
ternal coupling strength A = A' 3 ' = 0.01; bath temperatures 
T(B1)=0.4, T(B2)=0.2. 



any energy exchange (H^ R \ C NR = 1, H (F \ C F = 1). In 
this special case we find a non vanishing temperature pro- 
file independent of the initial state. In FIG.0]we show the 
temperature profile for a iV = 4 (N — 6) spin system for 
A = A^ B ' = 0.01. For this special coupling type, we observe 
the same profile even for larger A = A' B ) <= 10. In case 
of stronger internal coupling (A^ B -* = 0.01) we find again a 
linear profile but with a smaller gradient. Summarizing in 
all cases of Forster coupling together with a non-resonant 
coupling a non vanishing temperature gradient is reached, 
independently of the internal coupling strength. 

Finally, the internal interactions of the system are taken 
to be random next neighbor couplings. In FIG. |5] temper- 
ature profiles for different internal coupling strength A are 
shown. For decreasing A the profiles approach a linear de- 
pendence with finite gradient. 

Increasing A leads to a deviation from this linear pro- 
file: In case of very strong internal couplings the system is 
strongly entangled and therefore we find the spins locally 



in a totally mixed state. This behavior is similar for any 
random interaction. 

The distance between the expected stationary state 
and the actual final state of the system in dependence of 
internal coupling strength A can be found in FIG. = 
4). Again, the distance increases for a stronger internal 
coupling. But even for very small interactions we do not 
get exactly the expected stationary state: Any specific ran- 
dom coupling leads to a somewhat different but approxi- 
mate linear temperature profile. 



5 Flux 

In this section we examine the properties of the energy 
flux through the considered spin chain. We consider the 
change of energy of the single spin /i using the Liouvillc 
equation JJJ). 

(LjTrv {^} |1„) = (V|Tr„ {-i[H,p} + C™p} |1„>. 

Here, |1 M ) is the exited state of spin p and Tr^ means 
the trace over all other spins v [y ^ For a spin in the 
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Fig. 7. Fourier's Law: flux J over temperature gradient VT 
for a. N = 4 spin chain with Forster and non-resonant internal 
coupling (C*f = f , Cnr. = f , A = 0.1). 

J/10" 3 
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Fig. 8. The dependence of the flux J on the inverse chain 
length TV" 1 for C P = 1, Cnr = 1, A = 0.1 and N = 3, 4, 5, 6. 



middle of the chain (/j ^ 1 and /j, ^ N) the equation reads 
d(VlPlV) 



df 



2AC F (Im^-xiypHV-ilp) 

-MMWilflOpW)). (17) 



where trace over all the other spins is assumed. We inter- 
pret the quantity 



■W+i = 2AC F Im(l M M+ i|p|0 AI l At+ i) 



(18) 



as the flux from spin fi to spin fi+1 and the above equation 
shows that for the stationary state, the flux is the same 
for every /i, J = J MlM +i, as expected. 

We now test two important properties of the flux, 
its dependence on the temperature gradient and on the 
length of the chain. Again we consider now a system with 
a Forster and a non-resonant coupling. In FIG.EJwe show 
the flux J in dependence of the temperature gradient VT, 
the system fulfills Fourier's law 



const = J = kVT. 



(19) 



Only if VT is constant we can deduce the conductance k 
to be a constant material property over the homogeneous 
"wire". For the same T(B1), T(B2) the flux J is found to 
decrease with the chain length N: J cx (see FIG. |HJ), 
which underlines the fact that the conductance k is a bulk 
effect, not contact property. 



6 Conclusion 

We have studied a one-dimensional quantum system, a 
spin chain, coupled to heat baths within a Lindblad for- 
malism. For given local bath couplings, different interac- 
tions between the subsystems in the chain have been in- 
troduced, an energy transfer coupling (Forster coupling), 
a non-resonant coupling and an random "unbiased" inter- 
action. The chain properties have been characterized by 
the profile of the local mean energy and the energy-flux, 
J, respectively. 

In case of coupling to only one bath the system has 
been found to behave thermodynamical (i.e. flat temper- 
ature profile, controlled by bath) if we use the Forster 
coupling together with a non-resonant coupling. For ran- 
dom couplings we find a thermodynamical behavior for 
small interactions within the system only. 

The scenario for heat conduction, i.e. the spin chain 
coupled between a reservoir with high temperature and 
a reservoir with low temperature, shows different behav- 
ior in dependence of the type of the internal interactions. 
For energy transfer coupling only the heat conductivity 
diverges. With an additional non-resonant interaction we 
find a non vanishing temperature gradient. For longer 
chains the system shows the expected scaling properties. 
In case of small random interactions we always get a non 
vanishing temperature gradient, but for increasing interac- 
tion strength local thermodynamical behavior is violated. 

Even a class of small systems with just four spins may 
thus show local thermodynamical behavior under quan- 
tum mechanical modeling. For weakly coupled systems 
the emergence of a non vanishing temperature gradient 
appears to be generic if and only if the internal interaction 
contains some non-resonant part. For systems with strong 
internal interaction, however, the local mean energy profile 
typically shows significant deviations from linear behavior. 
In this case the definition of a local temperature becomes 
ambiguous. However, we hope to overcome this weak cou- 
pling limitation by introducing a coarse graining, i.e. by 
grouping a number of subsystems together such that these 
larger blocks exhibit effective weak coupling. Then a lo- 
cal energy profile on this effective reduced scale could be 
defined for which we expect to recover thermodynamical 
behavior. This would give an answer to the question on 
what local scale a temperature or local conductance could 
reasonably be defined. 

There are special interaction models (Forster coupling 
together with non-resonant coupling) , for which the profile 
does not depend on the interaction strength, a thermo- 
dynamical behavior can always be observed. For energy 
transfer only (Forster coupling) thermodynamical behav- 
ior can never be found. In this way it is possible to char- 
acterize the behavior of a large class of small systems ac- 
cording to their internal coupling. 

Other open problems relate to the dependence of the 
heat conductivity k on temperature. In the classical phonon 
theory of heat conduction in insulators there are some ex- 
pectations on this dependence based on experimental re- 
sults 15 ; . We hope to find a corresponding behavior in a 
fully quantum mechanical treatment. 
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